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PREFACE 


This  report  Is  part  of  RAND's  continuing  interest  in  the  geo¬ 
physical  effects  produced  by  nuclear  explosions.  It  is  a  continuation 
of  the  work  in  RM-5738  on  some  aspects  of  the  acoustic  effects  of  such 
explosions.  The  work  was  sponsored  by  the  United  States  Air  Force  and 
the  Advanced  Research  Projects  Agency. 
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SUMMARY 


The  pressure  pulse  generated  In  an  Isothermal  atmosphere  by  an 
energy  source  at  the  ground  Is  calculated  from  an  Integral  representa¬ 
tion  of  the  pressure  field  previously  derived  by  the  authors.  The 
shape  of  the  signal  Is  shown,  as  a  function  of  time,  at  several  dis¬ 
tances  from  the  source  for  a  fixed  altitude,  and  at  several  altitudes 
at  a  fixed  lateral  distance.  The  first  signal  to  arrive  at  any 
location  Is  a  high-frequency  acoustic  wave,  followed  by  a  low-frequency 
acoustic-gravity  wave.  The  onset  of  the  latter  is  marked  by  a  sharp 
front,  or  caustic.  At  any  Instant  (after  the  arrival  of  the  caustic), 
there  are  three  principal  frequency  components  at  any  location,  the 
lowest  of  which  becomes  dominant  as  time  progresses.  It  is  shown  how 
the  qualitative  features  of  the  flow,  as  well  as  the  exact  location 
of  the  caustic,  can  be  obtained  from  kinematic  theory. 
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I.  INTRODUCTION 


A  great  deal  of  attention  has  been  devoted  to  the  subject  of 
acoustic-gravity  waves  in  the  earth's  atmosphere.  Historically, 
interest  in  this  area  was  first  aroused  by  the  worldwide  detection  of 
atmospheric  disturbances  of  natural  origin.  This  interest  has  been 
renewed,  in  recent  years,  by  the  detection,  in  widely  spaced  regions, 
of  atmospheric  waves  generated  by  nuclear  explosions.  One  aspect  of 
such  waves,  and  the  one  which  underlies  this  study,  is  the  existence, 
at  ionospheric  heights,  of  outward  traveling  disturbances  produced  by 
explosions  near  the  ground.  For  a  summary  of  the  experimental  and 
theoretical  work  pertinent  to  this  particular  aspect  of  the  problem, 
the  reader  is  referred  to  Row  (1967),  where  an  extensive  bibliography 
is  also  to  be  found. 

Despite  its  well-known  shortcomings,  the  isothermal  atmosphere, 
because  of  its  tractabllity,  has  received  a  great  deal  of  theoretical 
attention.  It  has  been  remarked  often  in  the  literature  that,  al¬ 
though  caution  must  be  exercised  in  comparing  experimental  observa¬ 
tions  with  the  predictions  of  so  simple  a  model,  useful  insight  into 
the  propagation  of  atmospheric  disturbances  is  nevertheless  provided 
thereby.  One  of  the  principal  results  of  the  work  presented  here, 
viz.,  the  existence  of  a  front  (caustic)  bounding  the  spatial  domain 
of  the  low  frequency  part  of  the  disturbance,  could  be  considered  a 
case  in  point.  It  may  reasonably  be  expected  that  this  is  a  property 
which  the  real  atmosphere  shares  with  the  isothermal  atmosphere,  with 
the  shape  and  location  of  the  caustic  depending,  of  course,  on  the 
actual  density  variation  of  the  atmosphere. 
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Any  theoretical  treatment  of  the  atmospheric  waves  generated  by 
a  nuclear  explosion  obviously  requires  some  simplifying  assumptions 
about  the  source.  (For  a  discussion  of  this  aspect  of  the  problem, 
see  Pierce,  1968.)  In  this  paper,  the  interest  is  in  the  dispersive 
waves  far  from  the  source,  with  periods  much  larger  than  the  charac¬ 
teristic  period  of  the  source.  For  this  purpose,  an  Instantaneous 
point  source  seems  adequate.  In  particular,  a  point  energy  source  is 
assumed,  since  that  seems  more  appropriate  to  the  initial  conditions 
represented  by  a  nuclear  explosion  than  the  point  mass  and  point  im¬ 
pulse  sources  which  have  appeared  in  other  treatments.  Some  charac¬ 
teristics  of  acoustic -gravity  waves  depend  only  on  the  properties  of 
the  medium,  and  not  on  its  excitation.  However,  other  properties,  such 
as  the  spectrum  at  a  given  location,  obviously  do  depend  on  the  source. 
Such  properties  also  depend  on  the  boundary  conditions,  for  which 
reason  the  presence  of  the  ground  is  taken  explicitly  into  account  in 
this  work. 

In  the  following,  two  alternative  approaches  to  the  problem  dis¬ 
cussed  above  will  be  presented.  The  first  is  a  formulation  as  an 
initial  value  problem,  with  the  appropriate  boundary  conditions,  from 
which  the  characteristics  of  the  disturbance  can  be  calculated  by  the 
method  of  stationary  phase.  Details  of  this  method  have  been  presented 
elsewhere  (Cole  and  Grelfinger,  1968),  and  will  be  briefly  reviewed 
here  for  the  sake  of  completeness.  It  will  then  be  shown  how  the 
main  properties  of  the  dispersive  waves  in  the  asymptotic  field  can  be 
obtained  from  the  dispersion  relation  and  the  kinematic  formulation  of 
Whitham  (1961).  Thus,  in  addition  to  providing  useful  insight  into 
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che  propagation  of  atmospheric  disturbances,  the  isothermal  atmosphere 
also  furnishes  an  illuminating  example  of  the  relationship  between 
these  two  alternative  methods. 
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II.  INITIAL  VALUE  FORMULATION 


As  mentioned  In  the  Introduction,  this  section  Is  intended  prl- 
narily  as  a  brief  outline  of  the  treatment  by  the  authors  of  the  Initial 
value  problem,  the  details  of  which  appear  in  Cole  and  Greifinger 
(1968).  The  problem  considered  is  the  motion  in  an  isothermal  atmosphere, 
above  a  ground  plane,  produced  by  the  instantaneous  release  of  energy 
at  a  point  on  the  ground  (the  origin  of  Fig.  1).  For  an  isothermal 
atmosphere  in  a  uniform  gravitational  field,  the  equilibrium  conditions 
are  characterized  by  the  usual  exponential  distributions 


Vz)  p0(z) 

*  “  * 

P  P 


(1) 


where  P  and  {>  are  the  pressure  and  density  in  the  plane  z  -  0,  with 

the  positive  z  direction  taken  upward.  The  scale  height  h  is  given 

* 

in  terms  of  the  temperature  T  ,  the  specific  heat  ratio  y, and  the 

*  . 

gravitational  constant  g  by  the  relation  h  ■  RT  /g.  The  lsentroplc 

*  *  I/O  1/2  1/2 

sound  speed  is  c  ■  (yRT  )  ■  y  '  c  ,  where  c  »  (gh)  is  the 

8  8 

gravity  wave  speed. 

The  basic  parameter  of  the  problem  is 


e 


(y-DQc 

*  3 
P  hJ 


(2) 


where  Qq  is  the  energy  released  at  t  •  0.  This  parameter  roughly 
measures  the  energy  release  compared  to  the  internal  energy  in  a  scale 
height  volume.  For  sea  level  explosions,  where  the  scale  height  is  of 
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Fig.  1  -  Co-ordinate  system,  with  the  origin  at  the  point  of 
energy  release. 
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the  order  of  10  km,  e  ztz  1.5  x  10  ^  Y  ,  where  Y  is  the  hydrodynamic 

H  H 

yield  in  KT.  Thus,  e  «  1  for  explosions  of  up  to  several  MT  or  so. 
*  *  * 

With  h,  h/c  ,  P  and  p  as  the  units  of  length,  time,  pressure, 
and  density,  respectively,  an  expansion  in  the  small  parameter  c  is 
carried  out.  This  represents  the  flow  as  small  changes  superimposed 
on  the  ambient  state,  viz., 

P  ■  e  z  (1  +  cp+  . . . ) 

p  ■  e  2  (1  +  eo+  ...)  (3) 

•+  -+ 
q  ■  ev  +  ... 

V  ■  (u,w) 


where  u  is  the  radial  component  of  the  velocity  and  w  is  the  vertical 
component.  A  linearized  system  is  obtained  for  (p,  o,  v) ,  with  a 
6-function  energy  source. 

A  formal  Integral  representation  of  the  solution  of  the  equations 
of  notion  is  constructed  by  a  combination  of  Laplace  and  Hankel  trans¬ 
formation.  For  the  pressure  field,  for  example,  this  has  the  form 

00 

p(r,z,t)  -  e*^2  I  JQ(kr)  F(k;z,t)  k  dk  (4) 
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J  B 

2..  2 
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s  ds 


(5) 


y(s,k)  -  J  [ (s2+u2 (k) ) (s2+u2 (k) ) ] 1/2 


Re(y)  >  0 
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inll 


1/2 


[(i  +  k2  +  23  k)1/2  +  (i  +  k2-  23  k)1/2] 

■  Brunt-VHisHlH  frequency  (dimensionless) 


where  r  Is  the  cylindrical  radius  (Fig.  1). 

The  pressure  field  Is  the  main  quantity  of  Interest,  but  similar 

expressions  are  readily  obtained  for  the  other  variables.  The  usual 

z  /2 

amplification  factor  e  appears  In  these  expressions.  The  properties 

-1|2 

of  the  medium  are  contained  In  the  exponential  factor  e  ,  while 
the  boundary  conditions  and  the  properties  of  the  source  are  contained 
In  the  form  of  the  square  brackets.  The  Initial  conditions  are  satis¬ 
fied  by  taking  the  path  of  integration  In  (5)  to  the  right  of  all 
singularities  In  the  s-plane. 

A  closed  form  approximation  to  the  pressure  field,  valid  In  the 
main  wave  zone  for  large  r  and  t,  can  be  obtained  by  asymptotic  Inte¬ 
gration  of  (4)  and  (5).  The  procedure  Involves  the  transformation  of 
(5)  to  a  Fourier  Integral  over  those  frequencies  which  can  propagate 
for  a  given  (real)  k,  followed  by  two  successive  applications  of  the 
method  of  stationary  phase.  The  method  requires,  in  addition,  that 
the  asymptotic  approximation  for  the  Bessel  function  of  large  argument 
be  used. 

The  transformation  Is  accomplished  by  a  suitable  deformation  of 
the  path  of  Integration  In  the  s-plane.  The  Integrand  of  (5)  has 
branch  points  at  s  ■  +  iu^(k)  and  s  ■  +  iw^k),  poles  at  s  ■  +  Ik, 
and  an  essential  singularity  (due  to  the  p  In  the  exponent)  at  s  ■  0. 

If  the  s-plane  is  cut  In  the  manner  shown  In  Fig.  2  and  the  phase  of 


.  2  -  Contour  for  carrying  out  the  integration  in  the  s-plane, 
with  the  location  of  poles,  branch  points,  and  essential 
singularity  as  shown. 
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y  suitably  defined,  the  integrand  becomes  a  single  valued  function  of 
s  on  and  inside  the  contour.  By  Cauchy's  theorem,  then,  the  integral 
in  (5)  is  equal  to  the  sum  of  the  integrals  around  the  branch  lines 
plus  the  contribution  from  the  residues  at  the  poles.  The  essential 
singularity  is  outside  the  contour,  and  thus  does  not  contribute  a 
residue. 

The  contribution  from  the  poles,  which  arises  from  the  presence 
of  the  ground,  is  easily  calculated,  and  requires  no  approximation  of 
the  Bessel  function.  The  resulting  preseure  field  has  the  form 


PGR 


irrJJL  z 

i-  (2-y)  e  Y  _ t _ 

2"  Y  (t2-r2)3/2  * 


(6) 


This  "ground  wave"  is  a  non-oscillatory  cylindrical  wave,  which  exists 
only  behind  the  spherical  acoustic  front. 

The  dispersive  waves  arise  from  the  integrals  around  the  branch 
lines.  A  change  of  the  order  of  integration  makes  it  possible  to 
express  these  as  integrals  over  positive  frequencies  in  two  frequency 
bands.  The  pressure  field  in  the  low-frequency,  or  "acoustic -gravity," 
band  is  given  by 


PAG 


G(u,k)dkdu? 


(7) 


while  the  pressure  field  in  the  high-frequency,  or  "acoustic,"  band 
is  given  by 
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k  JQ(kr) f 
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-  (^1  -  ip)ellJZ] 


k^(w)  ■  a) 


2  1 
“  "4 


2  «2 
a)  -g 


1/2 


For  r  sufficiently  large  that  the  asymptotic  approximation  for 
the  Bessel  function  Is  valid,  the  Integral  over  k  may  be  carried  out 
by  the  method  of  stationary  phase.  The  resulting  expressions  for  the 
pressure  field  are 


PAG  *  6 


z/2  cos  » 
2x2R 


Im< 


J  iii  (B  -ml 


2  2W1//  2.1/2 

-b))(l/4-m  ) 


lcos$ 
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(^+  lp8)e1(“t+flR)  -  (^  -  ip8)ei(ut_nR>)da,: 


_  z/2  cos  ♦ 

PA  2 

2ir  R 


Re 


no 

I 

M2 


u)(ui2-B2)(u)2-1/4)1/2 


D(u,4>) 


■1 


(■^  +  ips)ei(u,t'nR)  +  -  ly8)ei(a)t+nR)j>  dm) 


(10) 
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where 


2  23  2  12  2  2 

D(u,$)  “  (to  ~3  T  cos  $  +  ~B  )iii  Bin  $ 


p  *  cos 
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(1/4-m2)(B2-m2) 


/  2  2  2  .  | 
(a)  -0  cob  $)  J 


1/2 
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(e2  -032) 
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This  representation  is  an  expansion  in  spherical  waves,  with  a 
spherical  phase  and  group  velocity  depending  on  the  pole  angle  It 
is  interesting  to  note  that  the  exponents  in  the  spherical  wave  repre¬ 
sentations  (9)  and  (10)  are  identical  with  those  in  the  exact  free- 
space  Green's  function  obtained  by  Pierce  (196?)  -and  Row  (1967). 

One  final  point  worthy  of  note  is  that  (9)  and  (10)  both  contain 
the  factor  cos  <p.  Thus,  in  this  approximation,  the  dispersive  wave 
train  vanishes  at  the  ground,  where  the  signal  consists  entirely  of 
the  non-osclllatory  cylindrical  wave  given  by  (6). 
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III.  PROPERTIES  OF  ASYMPTOTIC  PRESSURE  FIELD 
AND  NUMERICAL  RESULTS 


The  pressure  field  associated  with  the  dispersive  waves  has  been 
expressed  In  terms  of  Integrals  of  the  form 

/  f(u>,$)ei(ut  -  nR)  du  (11) 

over  the  two  bands  of  propagating  frequencies.  For  large  R  (spherical 
radius)  and  t,  the  exponentials  are  rapidly  oscillating  functions  of 
u,  and  the  Integrals  therefore  permit  a  second  application  of  the 
method  of  stationary  phase.  It  Is  easily  shown  that  stationary  points 
can  exist  only  for  the  exponentials  with  the  minus  sign.  Such  points 
must  be  solutions  of 

t  -  0  <12> 

and  are  therefore  functions  of  R/t  and  $,  only.  Some  Important  proper¬ 
ties  of  the  flow  field  can  be  Inferred  from  the  distribution  of 
stationary  points  In  (R/t,$)-space.  This  has  been  discussed  In  some 
detail  by  Cole  and  Grelflnger  (1968),  and  will  be  briefly  summarized 
at  this  point  to  provide  a  better  understanding  of  the  numerical 
results  to  follow. 

Obviously,  all  points  of  stationary  phase  must  lie  behind  the 
acoustic  front,  i.e.,  inside  a  circle  of  radius  R/t  ■  1  in  the 
(R/t, ♦)- plane.  In  the  acoustic  band,  1/2  £  w  <  »,  there  is  a  single 
stationary  point  associated  with  each  point  in  the  (R/t,$)-plane 


-19- 


behlnd  the  acoustic  front.  Along  any  ray  In  this  plane,  l.e.,  for 
any  fixed  value  of  $,  the  frequency  of  stationary  phase  decreases 
monotonlcally  from  to  -*■  «  at  the  acoustic  front  to  u  ■  1/2  at  the 
origin. 

In  the  acoustic-gravity  band,  B  cos  4>  <_  u  <_  S,  the  situation  Is 
quite  different.  For  any  given  angle  d> ,  there  are  two  stationary 
points  for  all  R/t  smaller  than  some  maximum  value.  At  this  maximum 
value  of  R/t,  the  two  stationary  points  coincide,  while  for  larger 
values  of  R/t  (at  the  given  4>) ,  there  are  no  stationary  points.  Along 
a  ray  in  the  (R/t,<J>)- plane,  the  two  frequencies  of  stationary  phase 
diverge  from  their  common  value  at  coincidence,  one  increasing  toward 
B  and  the  other  decreasing  toward  B  cos  (p  as  R/t  tends  toward  the 
origin. 

The  situation  for  the  acoustic-gravity  band  is  identical  to  that 
which  arises  In  the  case  of  Incompressible  flow  in  a  density  strati¬ 
fied  liquid  considered  by  Mowbray  and  Rarity  (1967).  As  in  that  case, 
the  locus  of  double  (coincident)  points  of  stationary  phase  defines 
a  front,  or  "caustic,"  representing  the  onset  of  the  disturbance. 
Between  the  caustic  and  the  acoustic  front,  the  acoustic-gravity  band 
Is  exponentially  small,  and  the  main  contribution  comes  from  the 
acoustic  band. 

The  location  of  the  caustic  can  be  calculated,  parametrically, 
from  Its  definition  as  the  locus  of  double  stationary  points.  The 
result  is  shown  In  Fig.  3,  where  the  caustic  has  been  plotted  in 
cylindrical  (r/t,  z/t)-space.  Also  plotted  (dashed  lines)  are  lines 
of  constant  ♦  “  (cot  -  R SI).  These  are  lines  along  which  the  phase 


3  -  Location  of  the  caustic  in  (r/t ,z/t)-space.  The  dashed  lines 
are  lines  along  which  the  phase  is  constant  at  any  instant. 
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ls  constant  at  a  given  t,  and  Intersect  the  caustic  In  a  cusp.  Crests 
(or  troughs),  as  particular  lines  of  constant  phase,  are  therefore 
representable  by  such  curves.  However,  since  lines  of  constant  phase 
are  lines  along  which  <J>  decreases  with  time,  crests  are  not  stationary 
In  this  diagram,  but  move  across  It  obliquely  from  left  to  right.  The 
number  of  crests  In  the  diagram  Increases  linearly  with  time,  new 
crests  being  continuously  created  at  the  origin. 

From  the  above  considerations,  the  frequency-time  history  of  the 
signal  at  any  location  can  be  described  qualitatively.  It  Is  to  be 
remembered  that  the  frequency  of  stationary  phase  Is  the  principal 
frequency  associated  with  the  group  of  waves  at  the  location  In  ques¬ 
tion.  Furthermore,  a  given  point  (R,$)  In  physical  space  corresponds, 
as  time  Increases,  to  points  approaching  the  origin  along  a  ray  In 
the  (R/t,$)-plane.  The  variation  of  the  frequency  of  stationary  phase 
along  such  a  ray  thus  provides  the  frequency-time  history  of  the 
signal  at  the  point  (R,<|>).  Clearly,  this  history  Is  qualitatively 
the  same  for  all  locations  at  the  same  angle  $,  which  all  map  Into  the 
same  ray  In  (R/t,4)>space.  Thus,  the  first  disturbance  to  arrive  at 
any  location  (R,<|>)  Is  the  acoustic  spherical  front,  containing  the 
very  high  frequencies.  The  principal  frequency  in  this  part  of  the 
disturbance  decreases  with  time,  tending  asymptotically  towards  u  “ 

1/2.  Some  time  after  the  passage  of  the  spherical  front,  the  caustic 
arrives,  carrying  the  frequency  of  the  double  stationary  point  associ¬ 
ated  with  the  angle  $.  This  part  of  the  signal  then  splits  into  two 
frequencies,  one  increasing  toward  8  and  the  other  decreasing  toward 
B  cos  4>  as  time  progresses.  Thus,  the  signal  at  any  (R,<J>)  consists, 
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after  long  times,  of  the  three  angular  frequencies  1/2,  6,  and 
3  cos  4>. 

Such  considerations,  of  course,  cannot  provide  the  shape  of  the 
signal,  which  depends  on  the  relative  amplitudes  and  phases  of  the 
contributing  frequencies.  This  can  be  determined  only  by  actual 
evaluation  of  (9)  and  (10) .  These  integrals  have  been  calculated  by 
the  method  of  stationary  phase,  and  the  results  are  presented  in  Figs. 
4,  5,  and  6.  The  numerical  results  are  based  on  a  specific  heat 
ratio  of  y  »  1.4  and  an  atmospheric  scale  height  of  h  »  8  km,  values 
which  are  more  appropriate  to  sea  level  than  to  ionospheric  altitudes. 
However,  no  attempt  is  being  made  here  to  correlate  these  results 
vlth  observation,  the  purpose  being  primarily  to  illustrate  the  dis¬ 
persion  of  the  pulse  with  distance  from  the  source. 

In  Fig.  4,  the  three  separate  frequency  components  are  shown  for 
the  signal  at  a  typical  location.  It  is  apparent  that,  after  the 
arrival  of  the  acoustic-gravity  components,  the  dominant  frequency  is 
the  lowest  frequency  component,  i.e.,  the  component  that  tends  toward 
u>  »  6  cos  <|>  as  time  increases,  it  is  this  part  of  the  signal  which 
has  been  calculated  approximately  by  Row  (1967)  for  an  unbounded  iso¬ 
thermal  atmosphere.  Figure  5  shows  how  the  pulse  spreads  out  with 
lateral  distance  from  the  source  at  a  fixed  altitude,  while  Fig.  6 
illustrates  the  change  in  shape  of  the  pulse  with  altitude  at  a  fixed 
lateral  distance  from  the  source.  Although  the  principal  contribution 
to  the  signal  is  from  the  lowest  frequency  component,  the  shape  of 
the  signal  is  noticeably  affected  by  the  presence  of  the  higher  fre¬ 
quency  components.  The  break  in  each  curve,  which  occurs  at  the 
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arrival  of  the  cauatic,  is  a  consequence  of  the  inapplicability  of 
the  8tatic.iary  phase  approximation  in  the  vicinity  of  the  caustic.  A 
more  accurate  approximation  is  the  neighborhood  of  the  caustic  would 
show  a  large,  but  continuous,  change  in  amplitude  with  passage  of  the 
caustic. 


t - 10  ■  EUpsed  time  (min) 


Fig.  4  -  Resolution  of  the  pressure  amplitude  at  a  typical  location 
into  i’-s  three  components.  The  elapsed  time  measures  the 
time  after  the  arrival  of  the  acoustic  front. 
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Fig.  5  -  Response  of  an  isothermal  atmosphere  to  an  energy  source  at  the 
ground,  for  fixed  altitude  and  varying  range.  The  elapsed  time  measures 
the  time  after  the  arrival  of  the  acoustic  front  at  the  given  location. 
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t-t0  •  Elapsed  tlma  (min > 

Fig.  6  -  Response  of  an  isothermal  atmosphere  to  an  energy  source  at  the 
ground,  for  fixed  range  and  varying  altitude.  The  elapsed  time  measures 
the  time  after  the  arrival  of  the  acoustic  front  at  the  given  location. 
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IV.  KINEMATIC  DESCRIPTION  OF  DISPERSIVE  WAVE  TRAIN  AND  CAUSTIC 

In  this  section,  it  will  be  shown  how  the  overall  geometric 
features  of  the  dispersive  wave  pattern  in  the  asymptotic  field  can  be 
obtained  from  the  dispersion  relation  and  the  kinematic  theory  of 
Whitham  (1961).  In  fact,  although  it  will  not  be  demonstrated  ex¬ 
plicitly  here,  inclusion  in  the  theory  of  energy  considerations 
(Whitham,  1961)  makes  it  possible  even  to  determine  the  relative 
amplitudes  of  the  various  frequency  components.  Thus,  only  the  precise 
phase  relationships  are  not  directly  obtainable  from  the  theory.  In 
every  other  respect,  the  results  obtained  by  kinematic  theory  are  iden¬ 
tical  with  those  obtained  by  the  method  of  stationary  phase.  This  is 
not  surprising  in  view  of  the  close  connection  between  the  assump¬ 
tions  underlying  the  two  methods  (Lighthill,  1965).  However,  the 

equivalence  is  of  more  than  academic  interest,  since  kinematic  theory 

* 

seems  more  readily  adaptable  to  the  case  of  a  slightly  non-uniform 
medium  than  is  the  method  of  stationary  phase. 

As  far  as  kinematic  theory  is  concerned,  all  of  the  dynamics  is 
contained  in  the  dispersion  relation 

u>  ■  u  (k)  (13) 

which  is  a  functional  relationship  between  frequency  and  wave  number 
k.  The  principal  result  of  the  kinematic  theory  is  that,  for  a  homo¬ 
geneous  medium  with  a  dispersion  relation  of  the  form  (13),  groups  of 
waves  of  given  wave  number  (and  therefore  also  of  given  frequency) 

* 

The  isothermal  atmosphere  is,  of  course,  not  a  uniform  medium. 
However,  it  can  be  made  to  "look"  uniform  by  a  suitable  choice  of 
dependent  variables. 
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propagate  with  constant  group  velocity  c  , 

Is 

the  dispersion  relation 


given  by  the  gradient  of 


c 

g 


tu(ic)  . 


(14) 


Thu 8,  if  all  the  waves  originate  at  the  origin,  a  given  group  travels 

radially  outward  from  the  origin  with  constant  speed  (provided  no 

reflections  take  place  at  a  boundary,  as  would  be  the  case  if  the 

source  were  above  the  ground).  Thus,  the  group  velocity  is  c  ■  f/t, 

8 

where  is  the  position  of  the  group  at  time  t.  These  relations  make 

it  possible  to  determine  the  location  in  ($,t)  of  groups  of  given 

wave  number  and  frequency.  This  procedure  will  now  be  carried  out 

In  some  detail  for  the  case  in  question. 

It  can  be  shown  from  Cole  and  Greiflnger  (1968)  (cf  Pierce, 

—  _2  /  2 

1963)  that  the  homogeneous  wave  equation  for  the  quantity  p  -  e  p, 
where  p  is  the  pressure  perturbation  defined  by  (3),  Is 


A  2  2  2  2  2 

rJL  3 _  ,  3__  lv  .  2.  3  ,,  - 

<  A  •  2  ^  2”a  ”  2  +  ®  ^  2+  2^p  °* 

at  at  a*  at  ax  ay 


(15) 


For  plane  wave  solutions  of  the  form  exp[i(k  x  +  ky  +  kz-  ut)], 

x  y  z 

the  resulting  dispersion  relation  can  be  written 


y  (k2  +  k2  +  j)  +  [  j  (k2  +  k2  +  h2  -  82  k2]1/2 
2  r  z  A  —  A  r  z  A'  rJ 


.2  i2  a.  i2 

k  ■  k  +  k 
r  x  y 


(16) 
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The  plus  sign  Identifies  the  high  frequency  acoustic  branch,  while  the 
minus  sign  is  associated  with  the  low  frequency  acoustic-gravity  branch. 
The  kinematics  apply  separately  to  each  branch,  leading  to  the  possi¬ 
bility  of  coexisting  families  of  waves. 

For  the  purposes  of  calculation,  it  is  more  convenient  to  write 
the  dispersion  relation  in  the  form 


2  .  1  2, 

2  (kz  +  4  ‘  “  } 

J  2  2 

<r  -u ) 


(17) 


treating  k  as  the  dependent  variable  and  u  and  k  as  the  independent 
r  z 

variables.  In  this  form,  the  acoustic  branch  corresponds  to  frequen- 

2  2  1 

cies  oj  kz  +  and  the  acoustic-gravity  branch  to  frequencies 
2  2 

gj  <_  6  .  In  terms  of  these  variables,  the  components  of  the  group 
velocity  take  the  form 


r_ 

t 


r3(kJ)1 


U(»  )  J 

'a(kj) 

K 

z 

z  _ 

h*2z)\ 

t 

’a(k^)' 

(a) 

3(a)2) 

(18) 


Furthermore,  since  the  groups  travel  radially  outward,  the  direction 
(polar  angle)  $  in  which  the  group  propagates  is  given  by  tan  q>  ■  r/z. 
From  these  relationships,  it  is  possible  to  express  the  components 
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of  the  wave  number  of  a  given  group  in  terms  of  the  frequency  and 
the  angle  at  which  the  group  propagates;  thus, 


k 

r 


2,1  2.  .2 

to  -  to  )  sin  $ 

2  2  2  2  2 
.  (8  -  u  )(u  -  6  cos  *  ) 


1/2 


(19) 


k 


z 


12  2  2 
(±  -w  )(B  -w  )cos 

2  2  2 

(to  -8  cos  $) 


1/2 


The  requirement  that  k^.  and  kz  be  real  serves  to  define  the 
allowed  frequency  bands  for  groups  of  waves  propagating  in  a  given 
direction.  These  bands  are  shown  in  Fig.  7,  where  k^  and  k^  are 
plotted  against  to  for  a  given  $.  The  frequency  range  of  the  acoustic 
band  is  1/2  <_  to  <  <*>,  while  for  the  acoustic-gravity  band  the  range  is 
8  cos  <(>  <_  <*>  _<  6.  The  cutoff  frequencies  for  the  propagation  of  groups 
in  a  given  direction  are  thus  identical  with  those  obtained  by  the 
method  of  stationary  phase  described  earlier. 

With  kf  and  kz  given  by  (19),  it  is  now  possible  to  express  the 

2 

components  of  the  group  velocity  (18)  in  terms  of  u>  and  4>.  The 
resultant  group  velocity  is 


,,1  2W  2  2W  2  2  2  . .  1/2 
R  [(^  -w  )(8  “oj  )(co  “6  cos  $)] 

8  "  t  "  2  2.2  .2.1  X  .  2  . 

uj[(8  -w  )  +8  (•£  -8  )  sin 


(20) 


This  result  agrees  exactly  with  that  derived  by  Pierce  (1963)  from 
the  free-space  Green's  function. 


Wave  number 
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/9cos  £ 


Fig.  7 


Vertical  and  radial  wave  number  (units  of  1/h)  as  a  function 
of  angular  frequency  (units  of  c*/h) . 
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In  Fig.  8,  c  Is  plotted  against  m  for  a  given  value  of 
© 

Equivalently,  this  is  a  plot  of  R/t  against  u,  where  (R,t)  is  the 
location  in  physical  space  (at  the  angle  $)  of  the  group  of  frequency 
u>.  In  the  acoustic  band,  the  group  velocity  increases  monotonlcally 
from  0  to  1  as  the  frequency  Increases  from  1/2  to  °o.  Thus,  there 
is  one,  and  only  one,  frequency  in  this  band  at  each  location  between 
the  origin  and  the  acoustic  front.  In  the  acoustic-gravity  band,  on 
the  other  hand,  the  group  velocity  is  zero  at  both  ends  and  has  a 
single  maximum  in  between.  Thus,  there  is  a  maximum  R/t  <  1  to  which 
waves  in  this  band  can  propagate,  or,  equivalently,  a  fastest  group 
of  waves.  For  all  smaller  values  of  R/t,  there  are  two  groups  of  waves 
from  the  acoustic-gravity  band  (in  addition  to  the  single  group  from 
the  acoustic  band).  These  results  are  Identical  with  those  obtained 
by  the  stationary  phase  approximation. 

The  occurrence  of  a  front,  or  caustic,  appears  in  the  kinematic 
theory  as  a  maximum  in  the  group  velocity  in  any  direction.  The  condi¬ 
tion  for  the  caustic  is  therefore 

3c  (w,$) 

-Sr— ' 0  <21> 

which  serves  to  define  the  frequency  wc($)  at  the  caustic  in  a  given 
direction.  If  this  is  combined  with  (18),  an  equation  for  the  caustic 
can  be  obtained.  The  result  again  agrees  exactly  with  the  stationary 
phase  calculation  of  Cole  and  Grelfingcr  (1968). 

It  can,  in  fact,  be  demonstrated  more  generally  that  the  equation 
for  a  caustic  obtained  by  kinematic  theory  is  Identical  with  that 


-28- 


obtained  by  the  method  of  stationary  phase.  Since  the  group  velocity 
is  in  the  radial  direction,  the  condition  that  the  group  velocity  be 
a  maximum  in  a  given  direction  is  equivalent  to  the  condition  that  the 
derivative  of  the  group  velocity  be  zero  in  the  direction  of  the  group 
velocity.  The  condition  for  a  caustic  can  thus  be  written  (in  k-space) 


c 

g 


2 

c 


(22) 


where  the  first  equality  follows  from  the  usual  vector  relationships 
and  ^  x  c  ”0.  In  terms  of  Cartesian  components,  this  becomes 


i 


3k^8kj  3k j 


0 


(i  “  1,2.3)  . 


(23) 


The  system  (23)  for  the  three  components  -tt—  of  the  group  velocity 

J 

has  non-trivial  solutions  only  if  tho  determinant  of  the  coefficients 
vanishes,  i.e.. 


3  ui 

3ki3kj 


0  . 


(24) 


This  condition  for  a  caustic  is  identical  with  that  deduced  from  very 
general  considerations  of  stationary  phase  (Lighthill,  1965). 

In  the  kinematic  theory,  there  is  no  acoustic-gravity  signal 
ahead  of  the  caustic.  As  mentioned  in  the  preceding  section,  in  the 
stationary  phase  approximation  of  the  exact  integral  representation, 
the  signal  ahead  of  the  fastest  group  is  exponentially  small.  More¬ 
over,  a  more  sophisticated  evaluation  of  such  integrals  (Lighthill, 
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1965)  provides  an  accurate  description  in  the  immediate  neighborhood 
of  the  front.  The  kinematic  theory,  however,  does  give  the  correct 
qualitative  picture  of  the  dispersive  wave  trains  behind  the  front  as 
well  as  the  exact  location  of  the  front  itself.  In  that  it  can  be 
easily  generalized  to  a  medium  with  slowly  varying  properties,  it 
should  prove  a  useful  method  in  treating  such  problems. 
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